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O j| Abstract: Cavity array metamaterials (CAMs), composed of optical 

.^ . microcavities in a lattice coupled via tight-binding interactions, represent 

a novel architecture for engineering metamaterials. Since the size of the 
CAMs' constituent elements are commensurate with the operating wave- 
length of the device, it cannot directly utilise classical transformation optics 
F^ . in the same way as traditional metamaterials. By directly transforming the 

internal geometry of the system, and locally tuning the permittivity between 

^ . cavities, we provide an alternative framework suitable for tight-binding im- 

iS I plementations of metamaterials. We develop a CAM-based cloak as the case 

^ . study. 
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1. Introduction 

Metamaterials, as an arrangement of artificial elements designed to achieve novel electromag- 
netic properties, represent a major advancement in material science. Typically constructed with 
conductive elements smaller than the operating wavelength, these structures have demonstrated 
novel electromagnetic behaviour, such as negative refraction fT\ and cloaking f2l. Photonic 
crystals (PhCs) may also mimic metamaterial behaviour [3rSi • Tight-binding models have re- 
cently been investigated as a novel metamaterial architecture ||6ll3 . A realization of this is the 
cavity array metamaterial (CAM) - an optical lattice composed of evanescently-coupled micro- 
cavities 1 6]. The quantum mechanical properties of the metamaterial are produced by coupling 
each cavity to single atoms that can be individually controlled. By tailoring the system's disper- 
sion properties, it can show effects generally found in classical metamaterials, such as negative 
refraction, evanescent wave amplification, and anomalous waveguiding properties, but with the 
distinct exception that it operates fundamentally at the quantum level. Tuning the atomic prop- 
erties opens up the possibility of reconfigurable metamaterials with quantum properties. 

Pivotal in the development of metamaterials was the advent of transformation optics (TO) as 
a methodology for designing novel electromagnetic devices, paving the way to the realization 
of a cloaking device f2l . Critically, because C AMs operate on a wavelength commensurate with 
the size of its constituent elements, they cannot directly utilise classical TO in the same way as 
traditional metamaterial systems. Here we develop a framework for TO using CAMs. 

Traditional TO works by applying a covariant co-ordinate transformation to Maxwell's equa- 
tions and then finding the required permittivity, e, and permeability, /i, for light to follow 
geodesies in this transformed co-ordinate system which in the original co-ordinate system 
would be curved paths Om. An attempt has been made to find a quantum equivalent that 
involves an invariant transformation of the Schrodinger equation, predicting cloaking of matter 
waves |10|. As a fundamental technique this work may find applications in systems such as 
ultracold rubidium atoms in optical lattices where the wavelength of the atom is many times 
greater than the lattice spacing so the system can be treated as an effective medium, although 
there are significant technical difficulties that need to be overcome first. 

PhCs as a metamaterial have the advantage that electromagnetic waves can propagate through 
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Fig. 1. (Color online) (a) Schematic of a 2D array of coupled optical cavities that guides 
electromagnetic waves to form a cloaked region, (b) The inter-cavity coupling k can be 
achieved through evanescent field overlap, (c) An infiltrated silicon-based PhC. The refrac- 
tive index of the liquid crystal infiltrated into air pores (blue) can be tuned with an applied 
electric field, providing a mechanism to dynamically control local intercavity couplings. 



them with relatively low loss. Classically, PhCs can effect TO by either varying the dielectric 
composition Pl or geometrically transforming the PhC [5|. Our approach is in the spirit of 
the latter, in that we physically deform the medium. However our system differs from Ref. [5 1 
because we are working with a quantum mechanical platform. Unlike Ref. 1 10] we work with 
photons and not matter waves. In particular, we will look at the propagation of a single photon. 
In Sec. [2] we discuss the guiding of light in CAMs to effect TO, using the cloaking of a 
circular region as a case study. We then calculate in Sec. |3]the required permittivity distribution 
in a coupled optical array to implement the circular cloak. 



2. Cloaking 

In CAMs, light propagation can be directly controlled by the lattice geometry and intersite 
hopping frequency. Therefore, unlike the continuous medium of classical metamaterials, co- 
ordinate transformations can be directly represented by changing the lattice configuration. The 
goal then is to keep the Hamiltonian invariant whilst the lattice configuration changes according 
to the required propagation response. Thus TO applied here to CAMs is a literal transformation 
of the medium space. 

To demonstrate TO in CAMs we consider an array of coupled optical cavities, the simplest 
of the CAMs, and show how it can be used as a cloak. This is schematically shown in Fig.[TJa). 
This choice of system, which has become increasingly accessible in experiments ifTTIl . suffices 
to demonstrate the effects of modifying geometry. The methods here are transferable to other 
more complex CAMs such as the Jaynes-Cummings-Hubbard (JCH) system . 

Consider a two-dimensional (2D) array of interacting identical cavities, described by the 
Hamiltonian (h= 1), 

J^ = coY^ajai - K* ^ ajcij, (1) 

(iJ) 

where aj (at) is the photonic creation (annihilation) operator at site /, T.{ij) is the sum over all 
nearest-neighbor cavities, CO is the cavity resonance frequency, and K is the inter-cavity coupling 
rate. We have used a second quantized description as it provides a convenient way to express 
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Fig. 2. (Color online) Dispersion surface £^ of a square coupled-cavity lattice system in the 
first Brillouin zone. The frequency of the contours are labeled with values of {E — co)/k. 
The white contour indicates the operating frequency used to simulate the effective point 
source in Fig.[3la). The black contour indicates the operating frequency used to simulate 
the beam source in Fig.[3lc). 



the system, and it allows for a direct translation to more complex quantum mechanical systems, 
such as the JCH model. A classical electromagnetic description of the problem however, would 
be equally valid. Note that as the system is linear, and we are working explicitly with single 
photons, our results apply equally to many photons, where Hong-Ou-Mandel effects (T2\ are 
ignored. 

For a uniform square lattice, Bloch's theorem can be applied to find an analytical expression 
for the dispersion relation [fT3l . 



E = (0 — 2K[cos{kxd) + cos{kyd)] 



(2) 



where k = kjX + kyy is the wave vector associated with the crystal momentum and d is the 
distance between neighboring sites. The dispersion surface of the system is shown depicted in 
Fig. [21 The behavior of the propagating modes is then governed by the group velocity. 



:V-^E = 2Kd[^ix\{kxd)x- 



-^m{kyd)y\. 



(3) 



By spatially modifying the dispersion properties, or equivalently the shape of the isoenergy 
contours, of the array through the control of inter-cavity detuning or coupling, it has been 
shown that the light field can be guided to demonstrate unusual behaviors such as negative 
refraction ||6] and superprism and collimating effects [7 1. In contrast, here we deform the array 
to affect the spatially varying dispersion properties so that the system acts as a circular cloak. 
We apply the standard coordinate transformation for such a cloak O, i.e., in the new polar 
coordinates (/,0'), 

, b — a , 

r' = —-r^a, f = 0, (4) 

h 

for the region r G (0, b), where a is the radius of the hidden region and b — ai^ the radial width 
of the cloak. The transformation directly relocates the relative positions of the cavities within 
this region, as shown in Fig. Ob). 
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Fig. 3. (Color online) (a) A quasi-point source forms a spherical wave. The spherical wave 
bends around the circular cloak region, rendering it effectively invisible, (b) The position 
of the cavity sites in the circular cloak are transformed according to Eq. ^ with a = 5 and 
^ = 10, making an annulus feature, (c) Cloaking of a normally incident Gaussian pulse, 
superimposed in time to form a continuous beam. The circular cloak region is formed with 
a = 50, Z? = 100. The x and y axis indicate site co-ordinates. 



To illustrate the cloaking property of the system, we solve the time evolution of the system, 
"^(t) = —iJ^\i/{t), in the one-excitation (single-photon) manifold, with two different initial 
conditions. We use isofrequency contours and group velocities given by Eq. (2) and (3) to 
identify operating frequencies. In Fig. Oa), we initialize the excitation in a superposition of 
all the /:-modes at an isofrequency contour near the band edge (white curve in Fig. [2]), as this 
is where it is most circular, to mimic a point- like source: we plot the probability distribution 
|i/A(^)p. In Fig.Oc), the excitation is initialized as a Gaussian pulse in a single ^-mode for a 
square contour (black bold curve in Fig. O to minimize dispersion, and plot the superimposed 
distribution Y,t I V^(0 P- ^^^ t)oth cases, the photon propagates around the hole and continues as 
if the hole was not there. Any detector outside the annulus region of deformation would not 
be able to distinguish between the original system and the cloaking system, thereby allowing 
anything inside the hole to be hidden from view. 

Unlike traditional transformation optics where a change in local £ and /i are required to effect 
the co-ordinate transformation, in CAMs the co-ordinate transformation is a direct transforma- 
tion of the physical configuration of the system, and the inter-site properties are tuned to ensure 
the invariance of the Hamiltonian (Eq.[T]). In other words, the inter-cavity coupling for different 
inter-cavity spacing needs to remain the same for all sites / and j (Kij = K, \/ij). In the next 
section we show how this can be achieved by tuning the intercavity permittivity. 



3. Permittivity detuning 

Consider a ID coupled cavity system with individual cavity width w, and permittivity profile 
described by step-functions such that £{x) = £a inside cavities and 8^ in-between cavities, as 
depicted in Fig.fTJb). Solving Maxwell's equation for a single cavity V^E^ = —£{x)co^/c^Eq 
yields the fundamental cavity eigenmode with, 

iFofx)!-/^"^'^^^'')' \x\<w/2 

where the normalization constant A satisfies J £a{x)\Ea{x)\^dx = 1. In a tight-binding model, 
the nearest-neighbour coupling is given by lfT4l . 



/w/2 
\Ea{x)\\Ea{x-d)\dx. (6) 

-w/2 



To keep J^ invariant from lattice deformation, we solve Eq. © for a constant inter-cavity 
coupling while allowing e^ to vary. 

Optical microcavities such as microdisks, pillars and toroids in general, and 2D PhC res- 
onators in particular, represent a versatile platform for realizing this scheme. PhC cavities are 
created by introducing a local inhomogeneity in a periodic dielectric lattice. By placing such 
defects in close proximity, they become evanescently coupled to form a coupled-cavity array. 
Alternatively, this coupling can be achieved through a line-defect, which has demonstrated 
resonant transmission across over 100 cavities in a thin silicon PhC membrane liTTI . In both 
cases, the effective inter-site permittivity (e^) is related to the refractive indices of the dielectric 
medium and hole infiltrations |15|, and also the holes' size, geometry and position |[T6l in the 
crystal regions between the cavities. 

That the index of liquid crystal infiltrated into air pores of PhC can be altered substantially 
using an applied electric field |15| [Fig.IlJc)], represents a means to dynamically tune the lo- 
cal inter-cavity coupling. Mechanisms to change the index of the background material have 
also been developed including local heating ifTTl . chemical etching ifTSl . electrical carrier in- 
jections |[T9l , near-field tip perturbation [20|, photosensitive material illumination ll2T1l . carbon 
dot deposition (221, fiber- taper probing [23], ion-beam irradiation [24J, and nanoparticle intru- 
sion 1251 . Controlling the atomic-cavity detuning in coupled atom-optical cavities has also been 
shown to play a similar role to local index modification (SI . 

As a specific case study, we consider an infiltrated silicon-based PhC [Fig.[T];d)] implemen- 
tation operating at A = 1 .5 /im with typical setup characteristics of e^ = 1 1 .7, e^ ^ 2.3 (prior to 
transformation), /c = 10^^ rad/s, and w = A/2. To achieve the above cloaking prescription with 
uniform CO and £«, we require the distribution of the inter- site permittivity to take up values 
shown in Fig. (H for the example in Fig. [3] The permittivity in-between neighboring sites is 
assumed to be constant. Each plot point in Fig. |4t)) represents an inter-site position, and is 
characterized by a radial distance and a permittivity value. Near the center of the cloak the 
deformation stretches the physical distance between nearest neighbors. To maintain the same 
hopping frequency, the intersite permittivity is decreased. Near the outer radius of the cloaking 
annulus, the sites are squashed closer together, and correspondingly the intersite permittivity 
here is greater than the baseline permittivity outside the cloaking region. 

4. Conclusion and outlook 

CAMs offers the potential of metamaterials with quantum properties. Here we have shown how 
transformation optics is implemented in the context of CAMs: a spatial transformation of the 
physical sites and careful tuning of the intersite coupling. As a case study, we showed how a 
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Fig. 4. (Color online) (a) Implementation of the cloaking device shown in Fig. [3] using a 
lattice distribution of inter- site permittivity e^. £/, values are denoted by the colored dots 
between the lattice sites, (b) Variation of the inter- site permittivity as a function of the 
distance from the cloaking center, in a quadrant (other quadrants are symmetrically equiva- 
lent). Near the centre, the intersite permittivity is lower than the baseline permittivity as the 
spacing between sites are stretched; whereas near the outer edge of the cloaking annulus 
the permittivity is higher as the sites are squashed closer together. The blue dots indicate 
the permittivity distribution along co-ordinate x = 25.5. 



single photon device can cloak a circular region. This work provides the framework for future 
investigation into quantum mechanical behaviour such as superposition and entanglement, in 
this system. These types of quantum mechancal behaviour could, for example, be achieved by 
coherently manipulating the states of atoms or atom-like systems coupled to the cavity sites. 
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